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SUMMARY 


A  simple  algorithm  has  been  developed  that  describes  the  vertical  structure  of 
the  visible  extinction  coefficient  and  the  liquid  water  content  (LWC)  for  low 
visibility/low  ceiling  conditions.  The  algorithm  Is  based  on  a 
parameterization  of  several  sets  of  data  describing  the  vertical  structure  of 
hazes,  fogs,  and  low  stratus  clouds.  The  algorithm  needs  as  input  values  the 
value  of  the  visible  extinction  coefficient  oe(0.55pm)  or  LWC,  and  the  cloud 
ceiling  height.  The  form  of  the  algorithm  is 


og  or  LWC  =  A  exp{B  exp(Cz)) 


(6a) 


where  the  coefficients  A  and  B  are  determined  by  boundary  values  and  the 
initial  value  at  the  surface,  C  determines  the  increase  of  o  or  LWC  with 
altitude,  and  z  is  the  altitude.  The  visibility  can  be  readily  determined 
from  the  visible  extinction  coefficient  or  vice  versa. 

When  the  fog  is  so  thick  at  the  surface  that  the  cloud  ceiling  height  can  not 
be  determined  (i.e.,  the  sky  is  obscured)  or  one  is  already  in  the  cloud,  the 
value  of  the  coefficient  C  is  directly  related  to  the  increase  in  LWC  as  one 
rises  in  altitude  at  the  wet  adiabatic  lapse  rate  and  to  the  change  in  droplet 
size  distribution  associated  with  convective  processes.  When  a  haze  or  fog  is 
present  and  the  cloud  ceiling  height  zQ  can  be  determined,  the  coefficient  C 
is  related  to  the  cloud  ceiling  height  via  the  equation 


C 


ln[ 


ln(E/A), 
ln(D/A)J  • 


(7) 


where  A,  D,  and  E  are  determined  by  the  boundary  conditions  and  initial  values 
at  the  surface,  as  explained  and  defined  in  the  text. 

Thus,  the  vertical  structure  in  visibility  or  LWC  can  be  determined  for  low 
visibility/low  stratus  cases  on  the  basis  of  two  inputs:  surface  visibility 
and  cloud  celling  height.  The  results  can  be  extended  to  other  wavelengths 
via  numerous  wavelength  scaling  laws  that  have  been  empirically  derived. 

The  utility  of  this  algorithm  is  that  it  relates  the  vertical  structure  to 
surface  observations  of  visibility  and  cloud  ceiling  height.  These  latter 
quantities  are  commonly  recorded  in  meteorological  observations  taken 
throughout  the  world.  Therefore,  frequencies  or  probabilities  of  occurrence 
can  be  determined  for  the  vertical  structure  of  low  visibility/low  ceiling 
conditions.  The  impact  of  these  conditions  on  sensor  and  system  performance 
can  now  be  more  readily  calculated. 
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INTRODUCTION 


The  development  of  precision-guided  munitions  and  sophisticated  electro- 
optical  sensors  has  placed  new  emphasis  on  the  term  visibility. 
Traditionally,  visibility  has  referred  to  visual  estimates  of  the  range  within 
which  certain  objects  were  discernable  against  the  horizon.  Often  the 

visibility  is  measured  and  recorded  automatically  by  several  types  of 
visibility  meters.  In  either  case,  the  emphasis  has  been  upon  the  horizontal 
visibility.  Increasingly,  United  States  and  North  Atlantic  Treaty 

Organization  (NATO)  military  forces  are  relying  on  new  surveillance  and  weapon 
systems  whose  sensors  must  function  over  a  slant  path,  where  variations  in  the 
vertical  as  well  as  the  horizontal  visibility  are  important. 

In  low  visibility  situations,  due  either  to  haze  or  fog,  a  growing  body  of 
observations1"5  are  showing  that  the  measured  visibility  at  the  surface  is  not 
representative  of  conditions  a  few  hundreds  of  meters,  or  even  tens  of  meters, 
above  the  surface.  Thus,  "slant  path  visibility"  can  be  significantly 

different  from  "horizontal  visibility."  In  a  significant  fraction  of  the 
cases  the  visibility  becomes  worse  as  the  height  above  the  surface 

increases.  These  cases  are  be  of  special  concern  in  this  report. 

Based  on  numerous  vertical  profiles  of  droplet  size  distributions,  an 

empirical  model  of  the  vertical  structure  of  hazes  and  fogs  has  been 
developed.®  This  earlier  work  is  now  extended  in  this  report  for  low 
visibility/low  stratus  conditions.  Two  quantities  of  physical  significance 
are  selected:  the  extinction  coefficient  at  a  wavelength  of  0.55pm,  which 
directly  relates  to  the  visibility  (see  appendix  A),  and  the  LWC,  which 
relates  to  many  microphysical  properties  of  fogs  and  which  scales  directly  to 


‘Hoihjelle,  D.  L. ,  et  al ,  Balloon-borne  Aerosol  Particle  Counter  Measurement 
Mode  in  Wintertime  at  Gratenwohr,  West  Germany ;  tCUM-DR-Zb-ii,  05  Army 
Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM,  1976. 

’Pinnick,  R.  G. ,  et  al ,  Vertical  Structure  in  Atmospheric  Fog  and  Haze  and  Its 
Effect  on  IR  Extinction,  ASL-)r-uUiO,  uS  Army  Atmospheric  Sciences  Laboratory, 
White  Sands  Missile  Range,  NM,  1978. 

’Loveland,  R.  B.,  et  al ,  Atmospheric  Characterization  Measurements  for 
Copperhead  Ground  Fog  Experiment,  Internal  Report,  05  Army  Atmospheric 
Sciences  Laboratory,  White  Sands  Missile  Range,  NM,  1978. 

“Lindberg,  J.  D. ,  et  al ,  Vertical  Distribution  of  Fog  and  Haze  Near  Greding, 
Germany,  During  February  and  March  i9Gu,  internal  Report,  us  Army  Atmospheric 
Sciences  Laboratory,  White  Sands  Missile  Range,  NM,  1980. 

’Lindberg,  J.  D. ,  Compiler,  Early  Wintertime  European  Fog  and  Haze:  Report  on 
Meppen  80,  ASL-TR-0108,  US  Army  Atmospnerlc  sciences  Laboratory,  White  sands 
Missile  Range,  NM,  1982. 

•Duncan,  L.  D. ,  J.  D.  Lindberg  and  R.  B.  Loveland,  An  Empirical  Model  of  the 
Vertical  Structure  of  German  Fogs,  ASL-TR-0071,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM,  1980. 
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the  extinction  at  many  infrared  wavelengths.’  The  algorithm  that  is  developed 
describes  the  vertical  structure  of  the  extinction  coefficient  and  LWC  as  a 
simple  function  of  two  variables:  the  initial  value  at  the  surface  and  the 
cloud  ceiling  height. 

BASIS  AND  DERIVATION 

Some  detailed  data  on  the  vertical  structure  of  fogs  and  hazes  have  been 
gathered  in  the  Federal  Republic  of  Germany  (FRG)  on  several  different 
occasions.1  J  *  5  Droplet  size  distributions  in  the  0.23um  to  47ym  range  have 
been  measured  from  a  balloon-borne  instrument,  thus  yielding  vertical 
profiles.  LWC  and  extinction  coefficients  at  desired  wavelengths  can  be 
calculated  from  these  measured  droplet  size  distributions. 

The  vertical  structure  of  these  profiles  has  been  examined  previously  by 

Duncan  et  al , *  who  characterized  the  vertical  structure  in  the  form 


y  =  a'x  +  b‘  , 


’Pinnick,  R.  G. ,  et  al ,  Relationships  Between  IR  Extinction,  Absorption  and 
Liquid  Water  Content  of  Fogs"j  ASL-TR-0U37,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM,  1979. 


‘Hoihjelle,  D.  L.,  et  al ,  Balloon-borne  Aerosol  Particle  Counter  Measurement 
Mode  in  Wintertime  at  Grafenwohr,  west  Germany,  LCUM-UK- /b-3,  US  Army 
Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM,  1976. 

‘Pinnick,  R.  G. ,  et  al ,  Vertical  Structure  in  Atmospheric  Fog  and  Haze  and  Its 
Effect  on  IR  Extinction,  ASL-TR-00iu,  US  Army  Atmospheric  Sciences  Laboratory, 
White  Sands  Missile  Range,  NM,  1978. 

"Lindberg,  J.  D. ,  et  al ,  Vertical  Distribution  of  Fog  and  Haze  Near  Greding, 
Germany,  During  February  and  March  19HU,  internal  Report,  us  Army  Atmospneric 
Sciences  Laboratory,  White  Sands  Missile  Range,  NM,  1980. 

•Lindberg,  J.  D. ,  Compiler,  Early  Wintertime  European  Fog  and  Haze:  Report  on 
Meppen  80,  ASL-TR-0108,  US  Army  Atmospheric  sciences  Laboratory,  wnite  Sands 
Missile  Range,  NM,  1982. 


•Duncan,  L.  D. ,  J.  D.  Lindberg  and  R.  B.  Loveland,  An  Empirical  Model  of  the 
Vertical  Structure  of  German  Fogs,  ASL-TR-0071,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM,  1980. 


LOG  D(Z+  20) 


LOG  D  (Z) 


Figure  la.  Relationship  between  the  extinction  coefficient  (0.55pm)  at 

altitudes  z  and  z  +  20  m.  The  vertical  lines  are  the  error  bars 
for  the  data  (after  Duncan  et  als). 


LOG  D  (Z) 

Figure  lb.  Relationship  between  LWC  at  altitudes  z  and  z  +  20.  The  vertical 
lines  are  the  error  bars  for  the  data  (after  Duncan  et  al‘). 
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where  x  =  logioD(z),  y  =  1ogi«D(z  +  20),  a'  and  b'  are  coefficients  that  were 
chosen  to  fit  the  data,  and  0 ( z )  is  the  value  of  the  desired  variable  (e.g. 
LWC)  at  the  altitude  z;  0(z  +  20)  is  then  the  value  of  this  variable  at  an 
altitude  of  z  +  20  m.  Thus,  one  can  work  stepwise  from  the  surface  up  through 
the  cloud  boundary  layer.  Figures  la  and  lb  show  the  fits  of  equation  (1)  to 
the  data.4 

Two  points  should  be  noted  about  figures  la  and  lb  before  proceeding.  First, 
each  figure  shows  that  two  lines  must  be  used  to  fit  the  data  set.  The  point 
of  intersection  of  these  two  lines  does  have  physical  significance,  as  will  be 
shown  in  the  next  section.  Second,  the  coefficients  of  the  equations  for  the 
two  lines  have  values  such  that  there  is  a  point  at  which  y  =  x  (i.e., 
logi0D(z  +  20)  =  log10D(z)).  This  means  that  there  is  an  upper  and  a  lower 
bound  to  the  variables  and  their  vertical  profiles,  or  in  other  words  there  is 
a  point  where  the  vertical  profile  of  o0  or  LWC  no  longer  increases  (or 
decreases)  with  altitude.  Again,  the  physical  interpretation  will  be 
discussed  in  the  next  section. 

Since  the  variable  y  in  equation  (1)  is  simply  the  variable  x(z)  at  an 
altitude  of  z  +  20  m,  the  equation  can  be  rewritten  as 

x(z  +  20)  -  x(z)  =  (a1  -l)x(z)  +  b‘  .  ,£) 

where  x(z)  has  been  subtracted  from  each  side.  Dividing  each  side  by  20  m 
(i.e.,  az)  and  recognizing  that  x(z  +  20)  -  x(z)  is  ax,  equation  (2)  becomes 

ax  _  (a1  -  1)  „  b' 

—  ^  x  •  (3) 


which  can  be  recast  as  a  simple  differential  equation  of  the  form 


+  ax  =  b 


This  has  a  general  solution  of  the  form 


,  b,  -az  b 

x  *  (c  -  >  +  * 


•Duncan,  L.  D. ,  J.  D.  Lindberg,  and  R.  B.  Loveland,  An  Empirical  Model  of  the 
Vertical  Structure  of  German  Fogs,  ASL-TR-0071,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM,  1980. 


12 


where  c  is  a  constant  of  integration  determined  by  the  boundary  conditions, 
and  a  and  b  can  be  found  from  figures  la  and  lb  and  equations  (1)  and  (3). 

Since  x  =  logi«D(z),  the  final  form  of  the  solution  of  the  extinction 
coefficient  (oe(0.55um))  or  LWC  is 


igor  LWC  =  A  exp[B  exp(Cz)] 


(6a) 


where 


A  = 


(6b) 


B  =  (In  10)  (c  -  ~) 

a 


(6c) 


where  In  10  =  2.3026...,  and  finally 


(6d) 


The  constant  of  integration  c  represents  the  logt»  of  the  initial  or  starting 
value  of  oe  or  LWC;  ^represents  logi,  of  a  bounding  value  for  ae  or  LWC. 
Thus,  the  coefficient  B  may  be  recast  in  the  form 


B  =  In  D  -  In  A  =  ln[j) 


(6e) 


where  D  is  the  initial  or  starting  value  of  the  variable  and  A  is  the  boundary 
value.  (Note  the  shift  to  natural  logarithms  and  exponentials.) 

Tables  la  and  lb  give  the  values  of  the  coefficients  in  equation  (6a)  for  the 
vertical  profiles  of  the  extinction  coefficient  and  LWC  as  parameterized  in 
figures  la  and  lb.  Note  that  there  are  two  equations  In  each  figure,  one  for 
each  line  segment,  and  therefore  two  sets  of  values  for  each  variable. 
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TABLE  la.  EXTINCTION  COEFFICIENT 


.1 

oe(0.55nm)  =  A  exp[B  exp(Cz)]  (km  ) 


Range  of 

Appl  icabll  ity 

0.398  to  7.08  km'* 
Haze/Fog 

7.08  to  92.1  km-i 

Thick  Fog/Cloud 

A 

0.398  km” 

92.1  km”1 

B 

In  (D/A) 

In  (D/A) 

C 

0.0125  m”‘* 

-0.014  m”' 

D 

Initial  value 

initial  value 

*C  also  may  be 
table  2  In  the 

calculated  as  an  explicit  function  of  cloud  celling  height;  see 
next  section. 

TABLE  lb.  LIQUID  WATER  CONTENT 

LWC  =  A  exp[B  exp(Cz)]  (g/m*) 

Range  of 

8.66(-5)*  to  2 .34 ( -2)  g/m1 

2 . 34 ( -2)  to  4.64(-l)  g/m1 

Appl icabil  f  ty 

Haze/Fog 

Thick  Fog/Cloud 

A 

8.66(-5)  g/m1 

4.64(-l)  g/m1 

B 

In  (D/A) 

_Hr* 

In  (D/A) 

_  1 

C 

0.008  m 

-0.015  m 

D 

Initial  value 

Initial  value 

*8.66(-5)  =  8.66  x  10”’ 


**C  also  may  be  calculated  as  an  explicit  function  of  cloud  ceiling  height; 
see  table  2  in  the  next  section. 


PHYSICAL  INTERPRETATION 

It  Is  now  worthwhile  to  examine  the  two  sets  of  coefficients  for  equation  (6a) 
and  to  determine  the  physical  significance  of  the  Intersection  of  the  two 
lines  shown  In  figures  la  and  lb.  For  this  purpose  It  Is  most  useful  to 
examine  the  algorithm  for  LWC.  The  two  lines  In  figure  lb  Intersect  at  a 
value  for  the  LWC  of  0.0234  g/m1.  Analyses  at  the  US  Amy  Atmospheric 
Sciences  Laboratory  (ASL)  of  numerous  vertical  profiles  have  shown  that  a 
cloud  boundary  can  be  reasonably  defined  for  a  value  of  LWC  between  0.02  and 
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0.03  g/m1.  Independent  analyses  of  numerous  cases  of  thick  fogs  and  clouds 
have  also  shown  that  0.03  g/m1  is  value  representative  of  the  boundary  or 
transition  between  fog  and  a  thick  fog/cloud.* 

The  point  of  Intersection  of  the  two  lines  in  figure  lb  thus  represents  a 
transition  between  two  regions  of  droplet  growth.  In  the  lower  region, 
representative  of  hazes  and  fogs,  the  growth  of  droplets  is  governed  more  by 
aerosol  hydration  below  and  at  supersaturation,  radiative  cooling  of  droplets, 
and  heat  transfer  to  the  ground.  In  the  upper  region,  representative  of  thick 
fogs  and  clouds,  the  droplet  growth  and  hence  tWC  increase  is  governed  more  by 
convective  motion.  The  transition  takes  place  where  LWC  reaches  a  value  in 
the  range  of  0.02  to  0.03  g/m5,  which  represents  either  a  fog  sufficiently 
dense  and  of  sufficient  vertical  depth  that  convective  cells  can  form,10  or 
the  lower  boundary  of  a  stratus  cloud.  The  precise  value  of  0.0234  g/m1  Is 
not  significant  in  itself,  but  only  represents  an  average  of  several  sets  of 
data. 

Figure  2  shows  the  comparisons  between  LWC  as  a  function  of  altitude,  as 
calculated  from  equation  (6a),  and  LWC  increase  due  to  air  rising  at  the  wet 


Figure  2.  Comparison  of  LWC  as  determined  from  equation  (6a)  and  the  Increase 
in  LWC  due  to  a  parcel  of  air  rising  at  the  moist  adiabatic  lapse 
rate. 


•Rogers,  C.  W.t  and  J.  T.  Hanley,  An  Algorithm  for  the  Increase  of  Liquid 


_ _  _ _ -iqui 

Water  Content  with  Height  in  Fog  and  Water  Hazes,  Caispan  Report  No.  67i i-H-1 , 
Buffalo.HY',  mW. -  - - 
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adiabatic  lapse  rate  (see  appendix  B).  Several  curves  are  shown  for  the 
latter  case,  illustrating  the  effects  of  temperature  and  lapse  rate.  The 
reasonable  agreement  between  the  two  types  o£  curves  shows  that  equation  (6a), 
with  the  second  set  of  coefficients  fron  table  lb,  is  an  average 

representation  of  the  increase  in  LWi  due  to  saturated  air  rising  at  the  wet 
adiabatic  lapse  rate.  Physically,  this  means  tha  when  a  fog  becomes 

sufficiently  dense,  overturning  motions  in  the  fog  began,’  10  so  that  further 
droplet  growth  and  LWC  increases  are  governed  by  convective  processes  and  the 
wet  adiabatic  lapse  rate.  The  presence  of  a  cloud  already  indicates 
convective  motions.  Deviations  from  the  wet  adiabatic  lapse  rate  may  be 
expected  near  the  earth's  surface  and  the  boundaries  of  the  cloud. 

The  value  of  the  coefficient  A  in  equation  (6a)  is  determined  from  setting  y  = 
x  in  figures  la  and  lb.  For  the  upper  region,  this  yields  a  value  of  0.464 
g/m,  which  is  in  good  agreement  with  the  fact  that  upper  bounds  to  LWC  in 
stratus  clouds  are  on  the  order  of  0.4  to  0.5  g/m.11  Of  course,  the  values 
can  be  smaller  for  thinner  clouds  that  are  not  undergoing  strong  radiative 
cooling  on  top.  Again,  the  precise  value  of  0.464  g/m5  only  represents  an 
average  from  several  sets  of  data.  For  the  lower  region,  the  boundary  value 
is  8.66  x  10"  g/m5.  There  are  two  interpretations  regarding  the  significance 

of  this  value.  On  the  one  hand,  this  corresponds  to  visibilities  on  the  order 
of  7  to  10  km  (light  haze),  which  is  representative  of  a  haze  aerosol  that  is 
beginning  to  grow  as  the  relative  humidity  rises  above  80  percent.12  On  the 
other  hand,  this  value  also  represents  the  lower  sensitivity  of  the 
instrumentation  used  to  measure  the  droplet  size  distribution.  In  actuality, 
both  effects  are  probably  present,  and  this  value  thus  represents  a  practical 
bound  for  the  condition  of  lowered  visibility.  Once  again,  it  is  an  average 
from  several  data  sets. 

In  an  analogous  manner,  the  point  of  intersection  of  the  two  lines  in  figure 
la  for  the  0 . 55un  extinction  coefficient  can  be  related  to  changes  in  droplet 
size  distributions.  These  changes  are  in  turn  due  to  the  change  in  physical 
processes  responsible  for  the  two  representations  for  LWC  in  figure  lb. 
Analysis  of  the  vertical  profiles  of  LWC  and  calculated  extinction 
coefficients  from  the  data  obtained  near  Meppen,  FRG,  in  19805  shows  that  the 


’Roach,  W.  T.,  et  al ,  "The  Physics  of  Radiation  Fog:  I  -  A  Field  Study," 
Quart  J  Roy  Meteorol  Soc,  102:313-333,  1976. 

“Pilie,  R.  J.,  et  al ,  "The  Life  Cycle  of  Valley  Fog.  Part  II:  Fog 
Microphysics,"  J  Appl  Meteorol ,  14:364-374,  1975. 

1 ‘Pruppacher,  H.  R.,  and  J.  0.  Klett,  Microphysics  of  Clouds  and 
Precipitation,  D.  Reidel  Pub  Co,  Dordrecht,  Holland,  1978. 

1  *Shettle,  E.  P. ,  and  R.  W.  Fenn,  Models  for  the  Aerosols  of  the  Lower 
Atmosphere  and  the  Effects  of  Humidity  Variations  on  Their  Optical  Properties, 
AFGL-TR- 79-0214,  Air  Force  Geophysics  Laboratory,  Hanscomb  AFB,  MA,  1979. 

’Lindberg,  J.  D.,  Compiler,  Early  Wintertime  European  Fog  and  Haze:  Report  on 
Meppen  80,  ASl-TR-0108,  US  Army  Atmospheric  sciences  Laboratory,  wnite  Sands 
MissDe  flange,  NM,  1982. 
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altitudes  are  the  same  for  the  points  where  ae  equals  7.06  km'1  and  LWC  equals 
0.0234  g/m1.  In  a  similar  manner,  the  upper  bound  of  92.1  km-1  for 

represents  the  limiting  value  of  transmission  in  a  stratus  cloud,  while  the 
value  of  0.398  km-1  represents  the  lower  limit  to  extinction,  or  the  limiting 
value  of  visibility  in  a  light  haze,  beyond  which  the  algorithm  should  not  be 
applied.  As  before,  these  boundary  values  are  averages  from  several  sets  of 
data. 

Figures  3a  and  3b  show  plots  of  the  vertical  profiles  of  the  extinction 
coefficient  and  LWC  using  equation  (6a)  and  tables  la  and  lb.  Different 
initial  values  have  been  selected  to  show  the  altitude  dependence  of  the 
various  profiles.  The  values  of  og  and  LWC  represented  by  the  intersection 
points  of  the  two  lines  in  figure6  la  and  lb  are  shown  by  dashed  vertical 
lines.  The  point  where  the  vertical  profile  intersects  the  dashed  vertical 
lines  is  the  nominal  lower  cloud  boundary;  the  altitude  is  indicated  in 
parentheses.  The  vertical  profile  up  through  the  cloud  can  be  obtained  by 
adding  the  profile  to  the  right  of  the  dashed  line  into  the  profile  to  the 
left  at  the  cloud  boundary,  and  then  extending  it  up  to  the  top  of  the 
c 1 oud . 


VISIBII.TY  (km) 

10.  I»  <>•> 


Figure  3a.  The  vertical  profile  of  the  0.55ym  extinction  coefficient  for 
various  Initial  values.  The  profiles  can  be  extended  to  the 
right  of  the  dashed  line  by  adding  the  profile  starting  at 
oe  =  7.08  km-1  and  starting  at  the  altitude  of  the  cloud  boundary. 
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Figure  3b.  The  vertical  profile  of  LWC  for  various  initial  values;  similar 
to  figure  3a. 

In  this  current  formulation,  the  initial  value  of  op  or  LWC  determines  the 
height  of  the  lower  cloud  boundary.  In  practice,  it  is  the  cloud  ceiling 
height  and  the  surface  value  of  visibility  that  are  the  observed  variables. 
(In  some  rare  instances  LWC  is  also  measured  at  the  surface.)  The 
deterministic  nature  of  equation  (6a)  comes  from  the  fact  that  the  coefficient 
C  is  single-valued.  It  has  been  shown  previously  that  for  thick  fog/cloud 
conditions,  this  value  is  related  to  the  increase  in  LWC  along  a  wet 
adiabat.  For  haze/fog  conditions,  the  value  of  C  given  in  tables  la  and  lb 
represents  an  average  of  several  different  haze  and  fog  cases.  Because  the 
intersection  points  in  figures  la  and  lb  can  be  related  to  the  cloud  boundary, 
it  is  also  possible  to  express  the  value  of  C  as  a  function  of  the  observed 
cloud  ceiling  height  z  : 


where  E  is  the  value  of  the  extinction  coefficient  or  LWC  at  the  upper  bound 
of  the  range  of  applicability  (i.e.,  the  values  of  oe  or  LWC  represented  by 
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the  dashed  vertical  lines  in  figures  3a  and  3b)  and  other  coefficients  are 
defined  as  before.  Table  2  gives  the  values  that  are  appropriate  in  each 
case. 

Figures  4a  and  4b  illustrate  the  cases  where  one  initial  value  for  oe  and  LWC 
has  been  picked,  but  several  different  ceiling  heights  have  been  specified. 
The  solid  line  gives  the  average  vertical  profile  from  tables  la  and  lb. 
Thus,  by  the  use  of  equations  (6a)  and  (7),  the  vertical  structure  for  low 
visibility/low-lying  stratus  conditions  can  be  specified  by  the  initial  value 
at  the  surface  and  the  cloud  ceiling  height. 


VISIHII.TY  (km) 

Id.  1.0  0.1 


CT(.  (0.5.">  pm)  (km  ' ) 


Figure  4a.  The  vertical  profile  of  the  0.55pm  extinction  coefficient  for 

various  cloud  ceiling  heights.  The  solid  line  shows  the  averaged 
profile  using  table  la. 
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Figure  4b.  The  vertical  profile  of  LWC  for  various  cloud  ceiling  heights. 
The  solid  line  shows  the  averaged  profile  using  table  lb. 


TABLE  2.  COEFFICIENT  C  AS  A  FUNCTION  OF  CEILING  HEIGHT 


C  = 


In  (E/A) 
In  (D/A) 


■] 


Range  of 

Appl  icability 

Extinction  Coefficient 

0.398  to  7.06  km" 

Liquid  Water  Content 
8.66(-5)  to  2 . 34( -2 )  g/mJ 

zc 

Cloud  ceiling  height  (mi 

Cloud  ceiling  height  (m) 

E 

7.08  km"1 

0.0234  g/m* 

A 

0.398  km'1 

8.66(-5)  g/m* 

0 

initial  value 

initial  value 
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COMPARISONS  WITH  DATA  AND  EXAMPLES  OF  USE 


The  empirical  representation8  of  the  vertical  structure  shown  in  figures  la 
and  lb  was  based  upon  the  data  gathered  in  Grafenwohr,  FRG. 1  1  A  validation 
of  the  model  would  be  to  compare  the  algorithm  based  on  these  data  with  an 
independent  set  of  data  obtained  in  Meppen,  FRG. s 

Figure  5  shows  the  comparison  between  the  algorithm  for  LWC  and  the  data 
obtained  from  the  balloon-borne  particle  spectrometer  at  Meppen.  In  this 
instance,  the  fog  was  so  thick  that  nc  ceiling  height  could  be  obtained,  and 
the  initial  value  of  LWC  was  greater  than  the  dividing  line  value  of  0.0234 
g/m’.  Therefore,  the  right-hand  column  of  table  lb  is  used  to  obtain  the 
value  of  the  constants  used  in  equations  (6a)  through  (6e) ;  an  initial  value 
of  0.04  g/m*  is  estimated  at  the  surface  from  the  data.  Substituting  these 
values  into  equation  (6a),  one  has 


LWC  =  0.464  exp[-2 .451  exp( -0 .015z) ] 


(6a-l ) 


The  comparison  with  the  data  in  figure  5  is  quite  good.  In  general,  the 
agreements  are  very  good  for  several  cases  that  were  tested,  with  the  values 
calculated  by  the  algorithm  being  within  +50  percent  of  the  data.  These 
accuracies  are  within  the  accuracy  limits  that  can  be  placed  on  the  data 
itsel f.* 


‘Duncan,  L.  D. ,  J.  D.  Lindberg  and  R.  B.  Loveland,  An  Empirical  Model  of  the 
Vertical  Structure  of  German  Fogs.  ASL-TR-0071,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM,  1980. 

‘Hoihjelle,  D.  L.,  et  al ,  Balloon-borne  Aerosol  Particle  Counter  Measurement 
Mode  in  Wintertime  at  Grafenwohr,  west  Germany!  ECoM-DR-76-3,  05  Army 

Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM,  1976. 

’Pinnick,  R.  G. ,  et  al ,  Vertical  Structure  in  Atmospheric  Fog  and  Haze  and  Its 
Effect  on  IR  Extinction,  ASL-TR-QUIU,  US  Army  Atmospheric  Sciences  Laboratory, 
White  Sands  Missile  Range,  NM,  1978. 

’Lindberg,  J.  D. ,  Compiler,  Early  Wintertime  European  Fog  and  Haze:  Report  on 
Meppen  80,  ASL-TT.-0108,  US  Army  Atmospheric  sciences  Laboratory,  white  Sands 
Missile  Range,  NM,  1982. 


*J.  D.  Lindberg,  private  communcation. 
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Figure  6  shows  the  comparison  between  the  algorithm  and  the  measured  LWC  for  a 
case  where  a  cloud  ceiling  height  can  be  estimated.  The  visibility  at  the 
surface  was  estimated  at  3  km,  with  thin  fog  (or  heavy  haze)  being  recorded  as 
an  obstruction  to  vision.  The  cloud  ceiling  height  was  estimated  from  the 
data  as  the  point  where  the  dividing  line  value  of  0.0234  g/m*  intersected  the 
LWC  profile,  z-  =  150  m.  (It  was  night  when  this  profile  was  taken,  and  no 
reliable  visual  estimate  of  the  cloud  ceiling  height  could  be  made.)  Using 
the  right-hand  column  from  table  2,  and  an  estimated  initial  value  for  LWC  at 
the  surface  of  0.0001  g/m’,  equation  (7)  for  C  yields 


r  -  1  i  rln(0. 0234/0. 0000866) ->  _  n  n?dd  - 
C  "  150  m  1[ln(0. 0001/0. 0000866  ]  "  0,0244  m 


(7-1) 


This  value  of  C  is  now  used  along  with  the  Initial  value  at  the  surface  in  the 
left-hand  column  of  table  lb  to  define  the  coefficients  in  equation  (6a) 


LWC  =  0.0000866  exp(0.l439  exp(0.0244z) ) 


(6a-2) 


which  defines  the  LWC  profile  from  the  surface  to  the  lower  cloud  boundary. 
The  LWC  profile  inside  the  cloud  is  now  found  by  using  the  right-hand  column 
of  table  lb  and  0.0234  g/m1  as  the  initial  value,  yielding 


LWC  =  0.464  exp[-2.987  exp(-0.015z‘ )] 


(6a-3) 


where  z'  is  now  the  altitude  above  the  cloud  base.  Again,  the  comparison  of 
the  calculated  LWC  profile  with  the  data  of  figure  6  Is  quite  reasonable; 
comparisons  of  several  profiles  of  this  type  with  similar  cases  from  the  ASL 
Meppen  80  tests  likewise  showed  good  agreement. 
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Figure  6.  The  comparison  between  the  algorithm  for  the  vertical  structure 
of  LWC  and  data  for  the  case  of  thin  fog  and  low  stratus  clouds 


Figure  shows  the  comparison  of  the  algorithm  for  the  visible  extinction 
coefficient  with  that  computed  from  the  balloon-borne  data  used  to  calculate 
the  LWC  profile  in  figure  6.  The  estimated  ceiling  height  is  found  by 
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Figure  7.  The  comparison  between  the  algorithm  for  the  vertical  structure 
of  the  0.55pm  extinction  coefficient  and  data  for  the  case  of 
thin  fog  and  low  stratus  clouds 

_  1 

extending  the  dividing  line  value  of  7.06  kr  upward  until  it  intersects  the 
data  profile.  As  before,  the  estimated  cloud  ceiling  height  is  150  m.  Thf 
value  of  the  coefficient  C  is  calculated  by  using  an  initial  value  of  0.4  km" 
and  the  coefficients  from  the  left-hand  column  of  table  2,  yielding 


150  m 


ln[ 


ln(7  .06/0.3981-1 
Into. 4/0. 398)  J 


=  0.0424  m-' 


(7-2) 


This  value  of  C,  along  with  the  initial  value  at  the  surface,  is  used  in  the 
left-hand  column  of  table  la  to  define  the  coefficients  in  equation  (6a) 


ae(0.55um)  =  92.1  exp[-2.566  exp(-0.014z' )] 


(6a-4) 


24 


where  z1  is  the  altitude  above  the  cloud  base.  Figure  7  shows  that  the 
comparison  between  the  data  and  the  algorithm  is  quite  good. 


In  each  of  figures  5  through  7  the  level  of  the  cloud  top  could  be  determined 
from  the  data.  In  practice,  the  top  boundary  height  of  a  cloud  deck  is 
usually  not  known,  nor  does  the  algorithm  itself  determine  where  the  top  of 
the  cloud  should  be.  Numerous  examples  of  vertical  profile  data  from  the  ASL 
Meppen  80  field  program  showed  that  low-lying  stratus  clouds  have  thicknesses 
mainly  in  the  range  of  100  to  300  m.  As  a  rather  rough  average,  a  nominal 
cloud  thickness  of  200  m  can  be  selected;  this  value  allows  the  LWC  to  grow  to 
0.4  g/m5,  which  is  a  reasonable  upper  limit.  Therefore,  when  the  cloud 
thickness  is  not  known,  a  default  value  of  200  m  is  suggested  for  use  in  the 
algorithm.  In  practice,  the  actual  cloud  thickness  is  not  a  crucial 
quantity.  The  attenuation  at  the  cloud  base  increases  so  rapidly  with 
increasing  altitude  that  after  the  first  SO  to  100  m  of  cloud  thickness, 
essentially  no  transmission  in  either  the  visible  or  infrared  occurs. 

DISCUSSION  AND  CONCLUSIONS 

The  preceding  sections  describe  a  simple  algorithm  that  has  been  developed  to 
describe  the  vertical  structure  of  the  visible  extinction  coefficient  and  LWC 
for  low  visibility/low  ceiling  conditions.  The  visible  extinction  coefficient 
has  been  selected  because  of  its  direct  relationship  to  visibility,  which  is  a 
commonly  observed  and  recorded  meteorological  quantity,  and  because  of  the 
numerous  existing  scaling  laws  that  relate  the  visible  extinction  coefficient 
to  the  desired  infrared  extinction  coefficient.  LWC  has  been  selected  because 
of  its  relation  to  many  fog  and  cloud  microphysical  properties  and  because  it 
directly  scales  to  the  extinction  or  absorption  coefficients  at  several 
infrared  wavelengths. 

The  algorithm  is  applicable  in  cases  of  lowered  visibility  (<  7  km)  and  low  or 
obscured  ceilings.  The  algorithm  has  two  sets  of  coefficients,  each  set 
applicable  to  a  different  regime.  The  first  regime  is  for  hazes  and  fogs, 
where  the  visibility  ranges  from  7  km  down  to  approximately  0.4  km,  and  a  low 
cloud  ceiling  is  present.  In  this  regime  the  algorithm  needs  as  input  values 
the  visible  extinction  coefficient  or  LWC,  and  the  cloud  ceiling  height.  The 
second  regime  is  for  thick  fogs  (visibility  <  0.4  km),  where  the  sky  is 
obscured,  or  when  one  is  already  in  the  cloud.  Here  the  input  parameter  is 
either  the  visible  extinction  coefficient  or  LWC;  appropriate  boundary  values 
have  been  determined  from  the  data.  Because  the  depth  of  thick  fog  and  the 
vertical  thickness  of  low-lying  stratus  clouds  are  not  easily  obtainable 
quantities,  a  default  value  of  200  m  is  suggested  when  this  information  is  not 
otherwise  available. 

The  algorithm  should  not  be  used  when  the  visibility  is  greater  than  7  km, 
when  there  is  no  cloud  ceiling,  or  when  the  cloud  ceiling  is  above  1  km. 
There  is  such  a  paucity  of  data  for  these  cases  than  an  extension  to  these 
regimes  is  not  currently  warranted.  As  more  information  on  the  vertical 
distributions  of  aerosols  in  the  planetary  boundary  layer  and  the  visibility 
vertical  structure  becomes  available,  these  difficulties  will  be  overcome. 

The  vertical  profiles  used  to  develop  this  algorithm  are  from  tethered  balloon 
flights,  all  of  which  are  for  altitudes  of  less  than  800  m  (above  the  launch 
site).  There  is  a  lack  of  data  to  determine  what  the  upper  bound  on  ceiling 
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heights  should  be.  It  Is  recommended  that  If  the  cloud  celling  height  is 
above  1000  m  (i.e.,  nominally  above  the  planetary  boundary  layer),  then  the 
algorithm  should  not  be  used. 

The  case  of  surface  haze  layers  and  shallow  radiation  fogs,  which  can 
significantly  reduce  horizontal  visibility,  is  being  studied  further  at  ASL. 
In  this  case,  a  cloud  ceiling  Is  usually  not  present,  and  the  visibility 
improves  rather  than  degrades  with  Increasing  altitude.  While  the  overall 
frequency  of  occurrence  of  this  type  of  vertical  structure  Is  probably  larger 
than  that  for  the  case  where  the  algorithm  is  applicable,  the  Impact  on 
certain  types  of  system  performance  Is  less,  because  the  transmittance  along 
the  air- to- surface  slant  path  Is  greater  than  along  the  equivalent  horizontal 
surface  path.  Current  research  and  data-gathering  efforts  are  being  carried 
out  at  ASL  to  develop  an  algorithm  for  this  type  of  vertical  structure. 

In  summary,  an  algorithm  of  simple  analytical  form  has  been  developed,  which 
describes  the  vertical  structure  for  low  visiblllty/low  ceiling  conditions. 
The  utility  of  this  algorithm  is  that  it  relates  the  vertical  structure  to 
surface  observations  of  visibility  and  cloud  ceiling  height.  These  latter 
quantities  are  commonly  recorded  in  meteorological  observations  taken 
throughout  the  world.  The  visible  extinction  coefficient  and  LWC  have  been 
used  as  parameters  for  the  algorithm.  The  visibility  itself  or  extinction 
coefficients  at  other  wavelengths  could  have  been  used  just  as  well.  Appendix 
C  provides  a  succinct  set  of  tables  of  coefficients  for  the  algorithm  for  the 
parameters  of  visibility,  0.55pm,  4.0um,  and  10.6pm  extinction  coefficients, 
and  LWC. 
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APPENDIX  A 


VISIBILITY  DEFINITION 

A  brief  statement  on  the  use  of  the  term  "visibility"  is  in  order.  Three 
terms  are  commonly  used:  visibility,  visual  range,  and  meteorological 

range.  The  latter  two  have  precise  definitions,1  although  all  three  are  often 
used  interchangeably.  The  visual  range,  Ry,  may  be  defined  as 


R  =  —  in  -  , 

V  a  e  * 

e 


( A- 1 ) 


where  a-  is  the  extinction  coefficient  in  the  visual  wavelength  band,  C  is  the 
inherent  contrast  of  the  target  against  the  background,  and  e  is  the  threshold 
contrast  of  the  observer.  The  meteorological  range,  Rm,  is  defined  as  above 
for  the  case  where  a  black  target  is  against  the  background,  so  that  C=l,  and 
where  the  visual  contrast  threshold  is  taken  as  e=0.Q2,  a  near-optimum  value 
for  daylight  conditions.  Thus, 


d  =  _I  ln  B  3J12 

m  a  in  0.02  o 
e  e 


(A-2) 


The  meteorological  range  is  often  taken  as  the  "visibility,"  but  it  should  be 
clear  that  this  is  an  optimum  visibility. 

In  practice,  the  visual  range  is  a  more  useful  quantity  because  it  allows  for 
a  target/background  contrast  of  less  than  unity  and/or  a  threshold  contrast  of 
more  than  2  percent.  Many  visibility  meters  are  calibrated  on  the  basis  of  an 
assuned  "observer"  threshold  contrast  of  5  percent  (or  a  combination  of  target 
contrast  and  perceptual  threshold  contrast  so  that  C/e  =  r-7^).  The  visual 
range  would  be 


R„  =  —  In  ttac 

V  oe  0.05 


(A-3) 


The  visual  range  defined  In  this  manner  is  a  slightly  more  conservative 
estimate  than  the  meteorological  range  (Rv  -  0.766  1^,) ,  but  it  correlates 
better  with  the  visibility  reported  In  meteorological  observations.1  The  term 


‘McCartney,  E.  J.,  Optics  of  the  Atmosphere,  John  Wiley  and  Sons,  Inc.,  New 
York,  1976. 
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"visibility"  in  the  text  Is  used  in  this  sense  of  visual  range, 
from  extinction  coefficient  to  visibility,  such  as  in  figures  3a, 
were  made  using  equation  (A-3). 


Conversions 
4a,  and  7, 


*Gordon,  J.  I.,  Daytime  Visibility:  A  Conceptual  Review,  AFGL-TR- 79-0257,  Air 
Force  Geophysics  Laboratory,  Hanscomb  AFB,  MA,  t97i). 


APPENDIX  B 


COMPUTATION  OF  LIQUID  WATER  CONTENTS  AND  LIQUID  WATER 
CONTENT  LAPSE  RATES  FOR  MOIST  ADIABATIC  PROCESSES 

LWC  and  its  lapse  rates  were  obtained  by  computing  the  change  in  the  water 
vapor  saturation  mixing  ratio  (aRs)  over  a  specified  height  interval  (10  m) 
and  then  converting  that  change  to  LWC  and  aLWC/aZ.  The  equations  used  were 


LWC  =  a(LWC)  +  LWCQ 


(B-8) 


ALWC/aZ  = 

where  R$  =  water  vapor  saturation  mixing  ratio 
Z  =  height  (meters) 

T  *  absolute  temperature  (degrees  Kelvin) 


(B-9) 


irw  =  moist  adiabatic  lapse  rate  (•/•  0.6*C/100  m) 
g  =  gravitational  acceleration  9.806  m  s-1 

R  =  gas  constant  for  air  2.8704  x  10“ 1  J/g  K 

es  =  saturation  vapor  pressure  for  water  (millibars) 

L  =  latent  heat  of  condensation  2.485  x  10*  J/g 

R*  =  universal  gas  constant  8.3144  J/K  mol 

mv  =  molecular  weight  of  water  vapor  18.016  g  mol"1 

pa  =  air  density  (1,275  g/m1  at  P0  =  1000  mbar  and  T0  =  0°C 

(  )  =  values  at  the  Initial  reference  level  (l.e.,  P  =  1000  mbar  and 
selected  T  ,  LWCQ  *  0.0234  g/m1),  or  values  at°the  previous 
reference  level. 
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APPENDIX  C 

TABULARIZED  ALGORITHM  COEFFICIENTS 

This  appendix  contains  the  tables  of  coefficients  for  the  vertical  structure 
algorithm.  The  derivations  for  the  0.55um  extinction  coefficient  and  the  LWC 
are  given  in  the  text.  The  tabularized  values  for  the  visibility  are  taken 
directly  from  the  0.55um  extinction  coefficient  via  the  relationship  (see  also 
appendix  A) 


Vis 


3.0 

o  (0.55ym) 
e 


(C-l) 


The  tabularized  values  for  the  4.0um  and  10.6um  extinction  coefficients  are 
derived  from  the  parametrized  data  in  an  earlier  report  by  Duncan  and 
associates.* 


'Duncan,  L.  D. ,  J.  D.  Llndberg  and  R.  B.  Loveland,  An  Empirical  Model  of  the 
Vertical  Structure  of  German  Fogs,  ASL-TR-0071,  US  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  W,  1980. 
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TABLE  C-l.  VISIBILITY  (Vis) 
Vis  =  A  expCB  exp(Cz)] 


Range  of 

Appl  icabil  ity 

7.5  to  0.42  km 

Haze/ Fog 

0.42  to  0.03  km 

Thick  Fog/Cloud 

A 

7.54  km 

0.0326  km 

B 

ln(D/A) 

ln(D/A) 

C 

1  ,  Pln(E/A), 

z  ln[Tn7D7*V 
c 

-0.014  m~‘ 

D 

visibility  at 
surface 

visibility  at  surface 
(or  at  cloud  boundary: 

0.425  km) 

E 

0.425  km 

(not  used) 

zc 

cloud  celling 
height*  (meters) 

(not  used) 1 

‘If  the  cloud  ceiling 

height  is  greater  than 

1  km,  the  algorithm  should 

probably  not  be  used. 

’If  the  height  of  the  fog  top  or  the  cloud  thickness  is  not  known,  assign  a 
default  value  of  200  m  and  use  this  portion  of  the  algorithm  only  over  that 
interval . 
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TA8LE  C-2.  0. 55ym  EXTINCTION  COEFFICIENT 


ae(0.55unO  =  A  exp[B  exp(Cz)] 


Range  of 

0.398  to  7.08  km' 

7.08  to  92.1  km" 

Appl  icabil  ity 

Haze/Fog 

Thick  Fog/Cloud 

A 

0.398  km'1 

92.1  km'1 

B 

ln(D/A) 

ln(D/A) 

C 

1  , _rln(E/A)-, 

zc  lnCln(D/A)] 

-0.014  m'1 

0 

initial  value  at 

initial  value  at  surface 

surface 

(or  at  cloud  boundary: 

7.08  km‘ 

E 

7.08  km-' 

(not  used) 

zc 

cloud  ceiling  height1 

(not  used) 1 

(meters) 

‘If  the  cloud  ceiling 
probaoly  not  be  used. 

height  is  greater  than 

1  km,  the  algorithm  should 

*If  the  height  of  the  fog  top  or  the  cloud  thickness  is  not  known,  assign  a 
default  value  of  200  m  and  use  this  portion  of  the  algorithm  only  over  that 
interval . 
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TABLE  C-3.  4.0ym  EXTINCTION  COEFFICIENT 


oe(4.0um)  =  A  exp[B  exp(Cz)] 


Range  of 

Appl Icabil  Ity 

3.16(-2)*  to  3.32  km’1 
Haze/ Fog 

3.32  to  92.1  km"* 

Thick  Fog/Cloud 

A 

3.161-2)  km-' 

92.1  km-* 

B 

ln( D/A ) 

ln( D/A) 

C 

1  ,  rln(E/A), 

z  ln[m7Ai] 
c 

-0.014  m'1 

D 

initial  value  at 
surface 

Initial  value  at  surface 
(or  at  c|oud  boundary: 

3.32  km' 

E 

Zc 

3.32  km"1 

cloud  ceiling  height1 
(meters) 

(not  used) 

(not  used) 1 

*3 - 16 ( -2)  =  3.16  x  10-* 


1  If  the  cloud  celling  height  is  greater  than  1  km,  the  algorithm  should 
probably  not  be  used. 

1  If  the  height  of  the  fog  top  or  the  cloud  thickness  Is  not  known,  assign  a 
default  value  of  200  m  and  use  this  portion  of  the  algorithm  only  over  that 
Interval . 
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TABLE  C-4.  10.6um  EXTINCTION  COEFFICIENT 

go(10.6pm)  =  A  exp[B  exp(Cz)] 


Range  of 

1.0(-2)*  to  1.67  km"1 

1.67  to  100  km'1 

Applicability 

Haze/ Fog 

Thick  Fog/Cloud 

A 

0.01  km'1 

100  km'* 

B 

ln( D/A) 

ln(D/A) 

C 

1  ,  r  1  n(  E/ A )  -, 

zc  ln[ln(D/A)] 

0.0125  m~‘ 

D 

initial  value  at 
surface 

initial  value  at  surface 
(or  at  cloud  boundary: 

1.67  km-1) 

E 

1.67  km"1 

(not  used) 

zc 

cloud  ceiling  height1 
(meters) 

(not  used) 1 

*1.0(-2)  =  1.0  x  10-1 

‘If  the  cloud  celling  height  is  greater  than  1  km,  the  algorithm  should 
probably  not  be  used. 


lIf  the  height  for  the  fog  top  or  the  cloud  thickness  is  not  known,  assign  a 
default  value  of  200  m  and  use  this  portion  of  the  algorithm  only  over  that 
i nterval . 


TABLE  C-5.  LIQUID  WATER  CONTENT 
LWC  =  A  exp[B  exp(Cz)] 


Range  of 

Appl  icabil  ity 

8.66(-5)*  to  2.34{-2)  g/m! 
Haze/ Fog 

?.34(-2)  to  4.64(-l)  g/m» 
Thick  Fog/Cloud 

A 

B 

8.66(-5)  g/ms 
ln(D/A) 

4.64(-l)  g/m* 

ln( D/A) 

C 

1  ,  rln(E/A), 

z  lnClnTD7AT] 
c 

-0.15  m~ 

D 

initial  value  at  surface 

initial  value  at  surface 
(or  at  cloud  boundary: 
2.34(-2)  g/m5) 

E 

2.34(-2)  g/m5 

(not  used) 

zc 

cloud  ceiling  height1 
(meters) 

(not  used) * 

*8.66(-5)  =  8.66  x  10“’ 


‘If  the  cloud  celling  height  is  greater  than  1  km,  the  algorithm  should 
probably  not  be  used. 

*  If  the  height  of  the  fog  top  or  the  cloud  thickness  is  not  known,  assign  a 
default  value  of  200  m  and  use  this  portion  of  the  algorithm  only  over  that 
i nterval . 
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ATTN:  ATZHCD-CS 
Fort  Gordon,  GA  30905 

Commander 

US  Army  Signal  Center  A  Fort  Gordon 

ATTN:  ATZHCD-Q 

Fort  Gordon,  GA  30905 

USAFETAC/DNE 

ATTN:  M-.  Charles  Glauber 
Scott  AFB,  IL  62225 

Commander 

Air  Weather  Service 

ATTN:  AWS/DNDP  (LTC  Kit  G.  Cottrell) 

Scott  AFB,  IL  62225 

Commander 

Air  Weather  Service 

ATTN:  AWS/DOOF  (MAJ  Robert  Wright) 

Scott  AFB,  IL  62225 

Commander 

US  Army  Combined  Arms  Center 
A  Ft.  Leavenworth 

ATTN:  ATZLCA-CAA-Q  (Mr.  H.  Kent  Pickett) 
Fort  Leavenworth,  KS  66027 

Commander 

US  Army  Combined  Arms  Center 
A  Ft.  Leavenworth 

ATTN:  ATZLCA-SAN  (Robert  DeKInder,  Jr.) 
Fort  Leavenworth,  KS  66027 

Commander 

US  Army  Combined  Arms  Center 
4  Ft.  Leavenworth 

ATTN:  ATZLCA-SAN  (Mr.  Kent  I.  Johnson) 
Fort  Leavenworth,  KS  66027 

Commander 

US  Army  Combined  Arms  Center 
4  Ft.  Leavenworth 

ATTN:  ATZLCA-WE  (LTC  Darrell  Holland) 
Fort  Leavenworth,  KS  66027 

President 

USAARENBD 

ATTN:  ATZK-AE-TA  (Dr.  Charles  R.  Leake) 
Fort  Knox,  KY  40121 
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Commander 

US  Army  Armor  Center  and  Fort  Knox 
ATTN:  ATZK-CD-MS 
Fort  Knox,  KY  40121 

Commander 

US  Army  Armor  Center  and  Fort  Knox 
ATTN:  ATZK-CD-SD 
Fort  Knox,  KY  40121 

Aerodyne  Research  Inc. 

ATTN:  Dr.  John  F.  Ebersole 
Crosby  Drive 
Bedford,  MA  01730 

Commander 

Air  Force  Geophysics  Laboratory 
ATTN:  OPA  (Dr.  Robert  W.  Fenn) 

Hanscom  AFB,  MA  01731 

Commander 

Air  Force  Geophysics  Laboratory 
ATTN:  OPI  (Dr.  Robert  A.  McClatchey) 
Hanscom  AFB,  MA  01731 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

ATTN:  Dr.  T.  J.  Goblick,  B-370 

P.0.  Box  73 

Lexington,  MA  02173 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

ATTN:  Dr.  Michael  Gruber 

P.0.  Box  73 

Lexington,  MA  02173 

Raytheon  Company 

Equipment  Division 

ATTN:  Dr.  Charles  M.  Sonnenschein 

430  Boston  Post  Road 

Wayland,  MA  01778 

Commander 

US  Army  Ballistic  Research  Laboratory/ 
ARRADCOM 

ATTN:  DRDAR-BLB  (Mr.  Richard  McGee) 
Aberdeen  Proving  Ground,  MD  21005 


Comnander/DI rector 
Chemical  Systems  Laboratory 
US  Army  Armament  Research 
&  Development  Command 
ATTN:  DRDAR-CLB-PS  (Dr.  Edward  Stuebing) 
Aberdeen  Proving  Ground,  MD  21010 

Conmander/Di rector 
Chemical  Systems  Laboratory 
US  Army  Armament  Research 
A  Development  Command 
ATTN:  DRDAR-CLB-PS  (Mr.  Joseph  Vervier) 
Aberdeen  Proving  Ground,  MD  21010 

Commander/Di  rector 
Chemical  Systems  Laboratory 
US  Army  Armament  Research 
A  Development  Command 
ATTN:  DRDAR-CLY-A  (Mr.  Ronnald  Pennsyle) 
Aberdeen  Proving  Ground,  MD  21010 

Commander 

US  Army  Ballistic  Research  Laboratory/ 
ARRADCOM 

ATTN:  DRDAR-TSB-S  (STINFO) 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Artny  Electronics  Research 
A  Development  Command 
ATTN:  DRDEL-CCM  (W.  H.  Pepper) 

Adel  phi,  MD  20783 

Commander 

US  Army  Electronics  Research 
A  Development  Command 
ATTN:  DRDFL-CG/DRDEL-DC/DRDEL-CS 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

US  Army  Electronics  Research 
A  Development  Command 
ATTN:  DRDEL-CT 
2800  Powder  Mill  Road 
Adel phi,  MD  20783 

Commander 

US  Army  Electronics  Research 
A  Development  Command 
ATTN:  DRDEL-PAO  (M.  Singleton) 

2800  Powder  Mill  Road 
Adel phi,  MD  20783 
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Project  Manager 
Smoke/Obscurants 
ATTN:  DRDPM-Sf* 

(Dr.  Anthony  Van  de  Wal ,  Jr.) 
Aberdeen  Proving  Ground,  MC  21005 

Project  Manager 
Smoke/Obscurants 

ATTN:  DRDPM-SMK-T  (Mr.  Sidney  Gerard) 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Amy  Test  &  Evaluation  Command 
ATTN:  DRS^E-AD-M  (Mr.  Warren  M.  Baity) 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Test  &  Evaluation  Command 
ATTN:  DRSTE-AD-M  (Dr.  Norman  E.  Pentz) 
Aberdeen  Proving  Ground,  MD  21005 

Di rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-AAM  (Mr.  William  Smith) 
Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  "  DRXSY-CS  (Mr.  Philip  H.  Beavers) 
Aberdeen  Proving  Ground,  MD  21005 

Di rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-GB  (Wilbur  L.  Warfield) 
Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-GP  (Mr.  Fred  Campbell) 

Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Material  Systems  Analysis  Activity 
ATTN:  DRXSY-GP  (H.  Stamper) 

Aberdeen  Proving  Grounds,  MD  21005 

DI rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-GS 

(Mr.  Michael  Starks/Mr.  Julian  Chernlck) 
Aberdeen  Proving  Ground,  MD  21005 


Di rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-J  (Mr  James  F.  O'Bryon) 
Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-LM  (Mr.  Robert  M.  Marchettl) 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

Harry  Diamond  Laboratories 
ATTN:  Dr.  William  W.  Carter 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

Harry  Diamond  Laboratories 
ATTN:  DELHD-R-CM  (Mr.  Robert  McCoskey) 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

Harry  Diamond  Laboratories 
ATTN:  DELHD-R-CM-NM  (Dr.  Robert  Humphrey) 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

Harry  Diamond  Laboratories 
ATTN:  DELHD-R-CM-NM  (Dr.  Z.  G.  Sztankay) 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

Harry  Diamond  Laboratories 
ATTN:  DELHD-R-CM-NM  (Dr.  Joseph  Nemarlch) 
2800  Powder  Mill  Road 
Adel  phi,  MD  20783 

Commander 

Air  Force  Systems  Command 

ATTN:  WER  (Mr.  Richard  F.  Plcanso) 

Andrews  AFB,  MD  20334 

Martin  Marietta  Laboratories 
ATTN:  Jar  Mo  Chen 
1450  South  Rolling  Road 
Baltimore,  W)  21227 
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Commander  Or.  A.  D.  Belmont 

US  Army  Concepts  Analysis  Agency  Research  Division 

ATTN:  CSCA-SMC  (Mr.  Hal  E.  Hock)  Control  Data  Corporation 

8120  Hoodmont  Avenue  P.0.  Box  1249 

Bethesda,  MD  20014  Minneapolis,  MN  55440 

Director  Director 

National  Security  Agency  US  Army  Enqr  Waterways  Experiment  Station 

ATTN:  R52/Dr.  Douglas  Woods  ATTN:  WESEN  (Mr.  James  Mason) 

Fort  George  G.  Meade,  MD  20755  P.0.  Box  631 

Vicksburg,  MS  39180 

Chief 

Intelligence  Materiel  Development  Dr.  Jerry  Davis 

A  Support  Office  Department  of  Marine,  Earth 


'JS  Army  Electronic  Warfare  Laboratory 
ATTN:  DELEW-I  (LTC  Kenneth  E.  Thomas) 
Fort  George  G.  Meade,  MD  20755 

The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
ATTN:  Dr.  Michael  J.  Lun 
John  Hopkins  Road 
Laurel  1  MD  20810 

Dr.  Stephen  T.  Hanley 
1720  Rhodesia  Avenue 
Oxon  Hill,  MD  20022 

Science  Applications  Inc. 

ATTN:  Mr.  G.  D.  Currie 
15  Research  Drive 
Ann  Arbor,  MI  48103 

Science  Applications  Inc. 

ATTN:  Dr.  Robert  E.  Turner 
15  Research  Drive 
Ann  Arbor,  MI  48103 

Commander 

US  Army  Tank-Automotive  Research 
&  Development  Command 
ATTN:  DRDTA-ZSC  (Mr.  Harry  Young) 
Warren,  MI  48090 

Commander 

US  Army  Tank  Automotive  Research 
A  Development  Command 
ATTN:  DRDTA-ZSC  (Mr.  Wallace  Mick,  Or 
Warren,  MI  48090 


and  Atmospheric  Sciences 
North  Carolina  State  University 
Raleigh,  NC  27650 

Commander 

US  Army  Research  Office 

ATTN:  DRXRO-GS  (Dr.  Leo  Alpert) 

P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 
Commander 

US  Army  Research  Office 
ATTN:  DRXRO-PP  (Brenda  Mann) 

P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 
Commander 

US  Army  Cold  Regions  Research 
A  Engineering  Laboratory 
ATTN:  CRREL-RD  (Dr.  K.  F.  Sterrett) 
Hanover,  NH  03755 

Commander/Di  rector 
US  Army  Cold  Regions  Research 
A  Engineering  Laboratory 
ATTN:  CRREL-RG  (Mr.  George  Altken) 
Hanover,  NH  03755 

Commander 

US  Army  Cold  Regions  Research 
A  Engineering  Laboratory 
ATTN:  CRREL-RG  (Mr.  Roger  H.  Berger) 
Hanover,  NH  03755 

Commander 

US  Army  Armament  Research 
A  Development  Command 
ATTN:  DRDAR-AC  (Mr.  James  Greenfield) 
Dover,  NJ  07801 
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Commander 

US  Army  Armament  Research 
A  Development  Command 
ATTN:  DRDAR-TSS  (Bldg  #59) 

Oover,  NJ  07801 

Commander 

US  Army  Armament  Research 
&  Development  Command 
ATTN:  DRCPM-CAWS-EI  (Mr.  Peteris  Jansons) 
Dover,  NJ  078C1 

Commander 

US  Army  Armament  Research 
&  Development  Command 
ATTN:  DRCPM-CAWS-EI  (Mr.  G.  H.  Waldron) 
Dover,  NJ  07801 

Deputy  Joint  Project  Manager 
for  Navy/USMC  SAL  GP 

ATTN:  DRCPM-CAWS-NV  (CPT  Joseph  Mi  cell) 
Oover,  NJ  07801 

Commander/Di  rector 

JS  Army  Combat  Surveillance  A  Target 
Acquisition  Laboratory 
ATTN:  DELCS-I  (Mr.  David  Longinotti) 

Fort  Monmouth,  NJ  07703 

Commander/Di  rector 

US  Army  Combat  Surveillance  A  Target 
Acquisition  Laboratory 
ATTN:  DELCS-PE  (Mr.  Ben  A.  Di  Campli) 
Fort  Monmouth,  NJ  07703 

Commander/Di  rector 

US  Army  Combat  Surveillance  A  Target 
Acquisition  Laboratory 
ATTN:  OELCS-R-S  (Mr.  Donald  L.  Folanl) 
Fort  Monmouth,  NJ  07703 

Director 

US  Army  Electronics  Technology  A 
Devices  Laboratory 
ATTN:  DELET-DD  (S.  Danko) 

Fort  Monmouth,  NJ  07703 

Project  Manager 
FIREFINDER/REMBASS 

ATTN:  DRCPM-FFR-TM  (Mr.  John  M.  Blal’o) 
Fort  Monmouth,  NJ  07703 


Commander 

US  Army  Electronics  Research 
A  Development  Command 
ATTN:  DRDEL-SA  (Dr.  Walter  S.  McAfee) 
Fort  Monmouth,  NJ  07703 

OLA,  2WS  (MAC) 

Holloman  AFB,  NM  88330 
Commander 

Air  Force  Weapons  Laboratory 
ATTN:  AFWL/WE  (MW  John  R.  El  rick) 
Kirtland.  AFB,  NM  87117 

DI  ’■ector 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-SL 

White  Sands  Missile  Range,  88002 
Director 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-SL  (Dolores  Angulano) 

White  Sands  Missile  Range,  NM  88002 

Director 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TDB  (Mr.  Louie  Dominguez) 
White  Sands  Missile  Range,  NM  88002 

Director 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TDB  (Mr.  William  J.  Leach) 
White  Sands  Missile  Range,  NM  88002 

Di rector 

USA  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-TGP  (Mr.  Roger  F.  Willis) 
White  Sands  Missile  Range,  NM  88002 

Di rector 

Office  of  Missile  Electronic  Warfare 
ATTN:  DELEW-M-STO  (Dr.  Steven  Kovel) 
White  Sands  Missile  Range,  NM  88002 

Office  of  the  Test  Director 
Joint  Services  EO  GW  CM  Test  Program 
ATTN:  DRXDE-TD  (Mr.  Weldon  Findley) 
White  Sands  Missile  Range,  NM  88002 

Commander 

US  Amy  White  Sands  Missile  Range 
ATTN:  STEWS-PT-AL  (Laurel  B.  Saunders) 
White  Sands  Missile  Range,  NM  68002 
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Commander 

US  Army  RAO  Coordinator 
US  Embassy  -  Bonn 
Box  165 

APO  New  York  09080 

Grumman  Aerospace  Corporation 
Research  Department  -  MS  A08-35 
ATTN:  John  E.  A.  Selby 
Bethpage,  NY  11714 

Rome  Air  Development  Center 
ATTN:  Documents  Library 
TSLD  (Bette  Smith) 

Grifflss  AFB,  NY  13441 

Dr.  Roberto  Vagi Io-Laurin 
Faculty  of  Arts  and  Science 
Dept,  of  Applied  Science 
26-36  Stuyvesant  Street 
Mew  York,  NY  10003 

Air  Force  Wright  Aeronautical  Laboratories/ 
Avionics  Laboratory 

ATTN:  AFWAL/AARI-3  (Mr.  Harold  Geltmacher) 
Wright-Patterson  AFB,  OH  45433 

Air  Force  Wright  Aeronautical  Laboratories/ 
Avionics  Laboratory 

ATTN:  AFWAL/AARI-3  (CPT  William  C.  Smith) 
Wright-Patterson  AFB,  OH  45433 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CF-R  (CPT  James  M.  Watson) 

Fort  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CD-MS 
Fort  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CF-R 
Fort  Sill,  OK  73503 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  NOAA  Liaison  Officer 

(CDR  Jeffrey  G.  Carlen) 

Fort  Sill,  OK  73503 


Commandant 

US  Army  Field  Artillery  School 
Morris  Swett  Library 
ATTN:  Reference  Librarian 
Fort  Sill,  OK  73503 

Conmander 

Naval  Air  Development  Center 
ATTN:  Code  301  (Mr.  George  F.  Eck) 
Warminster,  PA  18974 

The  University  of  Texas  at  El  Paso 
Electrical  Engineering  Department 
ATTN:  Dr.  Joseph  H.  Plerluissi 
El  Paso,  TX  79968 

Commandant 

US  Army  Air  Defense  School 

ATTN:  ATSA-CD-SC-A  (CPT  Charles  T.  Thorn) 

Fort  Bliss,  TX  79916 

Commander 

HQ,  TRADOC  Combined  Arms  Test  Activity 
ATTN:  ATCAT-OP-Q  (CPT  Henry  C.  Cobb,  Jr.) 
Fort  Hood,  TX  76544 

Commander 

HQ,  TRADOC  Combined  Arms  Test  Activity 
ATTN:  ATCAT-SCI  (Dr.  Darrell  W.  Collier) 
Fort  Hood,  TX  76544 

Commander 

US  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-DA-L 
Dugway ,  UT  84022 

Commander 

US  Army  Dugway  Proving  Ground 

ATTN:  STEDP-MT-DA-M  (Mr.  Paul  E.  Carlson) 

Dugway,  UT  84022 

Commander 

US  Army  Dugway  Proving  Ground 

ATTN:  STEDP-MT-DA-T  (Mr.  John  Trethewey) 

Dugway,  UT  84022 

Commander 

US  Amy  Dugway  Proving  Ground 

ATTN:  STEDP-MT-DA-T  (Mr.  William  Peterson) 

Dugway,  UT  84022 
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Defense  Documentation  Center 
ATTN :  DDC-TCA 
Cameron  Station  Bldg  5 
Alexandria,  VA  22314 
12 

Ballistic  Missile  Defense  Program  Office 
ATTN:  DACS-BMT  (Colonel  Harry  F.  Ennis) 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Defense  Technical  Information  Center 
ATTN:  DDA-2  (Mr,  James  E.  Shafer) 
Cameron  Station,  Bldg  5 
■Vexandria,  VA  22314 

Commander 

‘J S  Army  Materiel  Development 
&  Readiness  Command 

ATTN:  DRCBSI-EE  (Mr.  Albert  Giambalvo) 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Commander 

:JS  Army  Materiel  Development 
&  Readiness  Command 
ATTN,  DRCIDC  (Mr.  James  Bender) 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Defense  Advanced  Rsch  Projects  Agency 
ATTN:  Steve  Zakanyez 
1400  Wilson  31 vd 
Arlington,  VA  22209 

Defense  Advanced  Rsch  Projects  Agency 
ATTN:  Dr.  James  Tegnelia 
1400  Wilson  Blvd 
Arlington,  VA  22209 

Institute  for  Defense  Analyses 
ATTN:  Mr.  Lucien  M.  Biberman 
400  Army-Navy  Orive 
Arlington,  VA  22202 

Institute  for  Defense  Analyses 
ATTN:  Or.  Ernest  Bauer 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Institute  for  Defense  Analyses 
ATTN:  Dr.  Hans  G.  Wolf hard 
400  Army-Navy  Drive 
Arlington,  VA  22202 


System  Planning  Corporation 
ATTN:  Mr.  Daniel  criedman 
1500  Wilson  Boulevard 
Arlington,  VA  22209 

System  Planning  Corporation 
ATTN:  COL  Hank  Shelton 
1500  Wilson  Boulevard 
Arlington,  VA  22209 

US  Army  Intelligence  *  Security  Command 
ATTN:  Edwin  Speakman,  Scientific  Advisor 
Arlington  Hall  Station 
Arlington,  VA  22212 

Commander 

US  Army  Operational  Test 
&  Evaluation  Agency 
ATTN:  CSTE-ED  (Mr.  Floyd  I.  Hill) 

5600  Columbia  Pike 
Falls  Church,  VA  22041 

Commander  and  Director 
US  Army  Engineer  Topographic  Laboratories 
ATTN:  ETL-GS-A  (Mr.  Thomas  Neidringhaus) 
Fort  Bel  voir,  VA  22060 

Director 

US  Army  Night  Vision  & 

Electro-Optics  Laboratory 
ATTN:  DELNV-L  (Dr.  Rudolf  G.  Buser) 

Fort  Bel  voir,  VA  22060 

Director 

US  Army  Night  Vision  & 

Electro-Optics  Laboratory 
ATTN:  DELNV-L  (Dr.  Robert  S.  Rodhe) 

Fort  Bel  voir,  VA  22060 

Di rector 

US  Army  Niqht  Vision  S 
Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Mr.  Joseph  R.  Moulton) 
Fort  Bel  voir,  VA  22060 

Director 

US  Army  Night  Vision  & 

Electro-Optics  Laboratory 
ATTN:  DELNV-VI  (Luanne  P.  Obert) 

Fort  Bel  voir,  VA  22060 


Oi rector 

US  Army  Night  Vision 

X  Electro-Optics  Laboratory 
ATTN:  DELNV-V1  (Mr.  Thomas  W.  Cassidy) 
Fort  Bel  voir,  VA  20060 

Di rector 

'JS  Army  Night  Vision  X 
Electro-Optics  Laboratory 
ATTN:  OELNV-VI  (Mr.  Richard  J.  Bergemann) 
Fort  Relvoir,  VA  2206C 

Di rector 

US  Army  Night  Vision  X 
Electro-Optics  Laboratory 
ATTN:  OELNV-VI  (Dr.  James  A.  Batches) 

Fort  Bel  voir,  VA  22060 

Commander 

'JS  Army  Training  X  Doctrine  Command 

ATTN:  ATCD-AN 

Fort  Monroe,  VA  23651 

Commander 

US  Amy  Training  X  Doctrine  Command 

ATTN:  ATCD-AN-M 

Fort  Monroe,  VA  23651 

Commander 

US  Army  Training  X  Doctrine  Command 
ATTN:  ATCD-F-A  (Mr.  Chris  O'Connor,  Jr.) 
Fort  Monroe,  VA  23651 

Commander 

IJS  Army  Training  X  Doctrine  Command 
ATTN:  ATCD-IE-R  (Mr.  David  M.  Ingram) 
Fort  Monroe,  VA  23651 

Commander 

US  Army  Training  X  Doctrine  Command 
ATTN:  ATCD-M-I/ATCD-M-A 
Fort  Monroe,  VA  23651 

Commander 

IJS  Army  Training  X  Doctrine  Command 
ATTN:  ATDOC-TA  (Dr.  Marvin  P.  Pastel) 
Fort  Monroe,  VA  23651 

Department  of  the  Air  Force 
OL-I,  AWS 

Fort  Monroe,  VA  23651 


Department  of  the  Air  Force 
HQS  5  Weather  Wing  (MAC) 

ATTN:  5  WW/DN 

Langley  Air  Force  Base,  VA  23655 
Commander 

US  Army  INSCOM/Quest  Research  Corporation 
ATTN:  Mr.  Donald  Wilmot 
6845  Elm  Street,  Suite  407 
McLean,  VA  22101 

General  Research  Corporation 
ATTN:  Dr.  Ralph  Zirkind 
7655  Old  Springhouse  Road 
McLean,  VA  22102 

Science  Applications,  Inc. 

8400  Westpark  Drive 
ATTN:  Dr.  John  E.  Cockayne 
McLean,  VA  22102 

US  Army  Nuclear  X  Chemical  Agency 
ATTN:  MONA-WE  (Dr.  John  A.  Berberet) 

7500  Backlick  Road,  Bldg  2073 
Springfield,  VA  22150 

Di rector 

US  Army  Signals  Warfare  Laboratory 
ATTN:  DELSW-EA  (Mr.  Douglas  Harkleroad) 
Vint  Hill  Farms  Station 
Warrenton,  VA  22186 

Di rector 

US  Army  Signals  Warfare  Laboratory 
ATTN:  DELSW-OS  (Dr.  Royal  H.  Burkhardt) 
Vint  Hill  Farms  Station 
Warrenton,  VA  22186 

Commander 
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